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The F18 HARV is a modified F18 aircraft which is capable of flying in the post-stall regime in order to 
achieve superagility. The onset of aerodynamic stall, and continued into the post-stall region, is characterized 
by nonlinearities in the aerodynamic coefficients. These aerodynamic coefficients are not expressed as analytic 
functions, but rather in the form of tabular data. The nonlinearities in the aerodynamic coefficients yield a 
nonlinear model of the aircraft’s dynamics. Nonlinear system theory has made many advances, but this area 
is not sufficiently developed to allow its application to this problem, as many of the theorems are existance 
theorems and that the systems are composed of analytic functions. Thus, the feedback matrices and the 
state estimators are obtained from linear system theory techniques. It is important, in order to obtain the 
correct feedback matrices and state estimators, that the linear description of the nonlinear flight dynamics 
be as accurate as possible. 

The nonlinear simulation is run under the Advanced Continuous Simulation Language (ACSL). The 
ACSL simulation uses FORTRAN subroutines to interface to the look-up tables for the aerodynamic data. 
ACSL has commands to form the linear representation for the system. This is a two step process. The 
first step is to trim the system, which is identical to trimming the aircraft. This involves calculating the 
trim input, u*, for a given trim state, x*, such that the derivatives* of the state vector, x, are zero. The 
state space for the simulation is x = (u, w, q t 0)^ The second st€p is to calculate the Jacobians of the state 
transisiton map, f(x) with respect to xat x*, and the input mi^ g(x,u) with respect to u at x* and u*. 
For analytic functions this is accomplished by calculating the partial derivatives of the functions, but since 
analytic functions do not exist, and ACSL can not perform symbolic calculations, the Jacobians are formed 
by numeric differentiation. ACSL uses the central difference method to perform the numeric differentiation. 

dfi(x) | = /t ( x* + hjQj) — fc (x* — 6j Gj ) 

dxj lx=x * 26 j W 

Where the perturbation vector is defined as 6 = (Ucos(a)/ 100, Usin(a)/ 100, 5%, 2°) T . This is used to 
calculate the A matrix for the linear representation, x = Ax + Bu with x = x - x* and u = u - u*. 
Similarly the B matrix is calculated by the numerical differentiation of g( ). The result of this is that the A 
and B matrices are not only functions of x* and u*, but also of 6 . 

The state transition map, f(x), and the input map g(x,u), are derived from the equations of motion 
for the system. A predominant part of these equations are based on the aerodynamic forces and moments. 
These forces and moments are calculated via the standard methods by using a nondimensional coefficient. 
These coefficients are the results of wind tunnel testing and they are stored in data files. The data points 
are generally functions of two variables, a and Mach number. Occasionally, a third variable is added such 
as control surface deflection. For this work, emphasis was placed on the a-Mach space variables. There are 
different sets of a vaules and Mach number values at which the coefficients are known. For example, C mo 
are known for a E {-12., -4., -2., 0., . . .} and Mach number E {0.2, 0.6, 0.8, 0.85, . . .} but C mq are known for 
a E {— 4.,0.,4., . . .} and Mach number E {0.2, 0.6, 0.8, 0.9, . . .}. To calculate the particular coefficient, with 
a and Mach number either equal or not equal to values in the respective sets, a second order Lagrange 
interpolation is performed. If the desired a and M are both equal to elements in the respective sets, then 
the Lagrange interpolation yields the coefficient from the table. This is accomplished first by bracketting 
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the a and Mach number, M , between availible t i t and Mj , such that, a, < a < o; +1 and A/ ; - < M < Mj+ x . 
Letting Ao — Oj+i - a* and A A/ = A/j + 1 — A/j, and be the coefficient C at, a* and A/j, the following 
expression yeilds (7(o, A/). 


C(a, A/) 




A / ) T f f / ,j -f. i ( i — a )( A/ — Mj ) 
AoA A! + 


— c»f)( A/j-n - M) + ('i+i t j + ](<x - otj)(M - Mj) 

AoA A/ 

It can be seen from the form of Kq. 2 that for a constant o or M , the plot for the resulting one variable function 
C( ) is a series of connected straight line segments. (Calculation of the various aerodynamic coefficients, not 
including coefficients in nr, and summation into the total force and moment calculation are perforin by the 
FORTRAN subroutine SFAFRRF. SFAFRRF -takes into account the various sets of a and M that are 
availible. 

When the system is linearized about some trim point and the eigenvalues are compared to those of 
the DMS simulation, some differences arise. The eigenvalues also exhibited a dependence on 6. With the 
need for accurate linear representations for the system, the following work addressed the problem of how to 
minimize the differences of the two systems. 

Due to the dependence on the perturbation size, <S, in the numerical calculation of the Jacobian, elim- 
ination of the numerical differentiation would address this problem and possibly assist in minimizing the 
differences between the two models. 'The concept was to express the tabular data in a known functional 
form, that could have its derivative calculated analytically by a FORTRAN subroutine. Various possible 
functional forms were examined, such as the fourth order Lagrange, full order Lagrange, and quadratic, 
cubic, quartic splines. These functions were C x or greater over the entire a-Mach space, being at least C 1 on 
the boundaries of different a-Mach cells. The full order Lagrange interpolation and the various splines pro- 
duced oscillations that were not desirable. The fourt h order Lagrange interpolation required a data window, 
which developed discontinuities in the derivatives at the window boundaries. Applying these techniques to 
the many data tables during each iteration would he computer intensive, and discontinuities would still exist 
in the derivatives on the cell boundaries. 

An intermediate approach has been followed. Numerical differentiation has been retained with a larger 
perturbation step size, and the data space for “sparcoly ’ distributed data has been augmented. The larger 
perturbation step size is opposite to the expected limit in calculating a derivative. The use of the increased 
perturbation size has the effect of decreasing the magnitude of the second derivative of the particular coeffi- 
cient. The perturbation slop size lias been set equal to one half of the minimum parameter step size, in this 
study 6 a = 1°. The a perturbation which is useful for operations with the aerodynamic data space, but a is 
not a state space variable but it is related to the state variable w. A function 8 W — f(6 a ) was calculated and 
included in the formulation ol the linear model. The augmenting of the data space also provides the for the 
minimization of the second derivative. The data space is augmented by calculating intermediate values for 
the particular coefficient by using a fourth order Lagrange interpolation. For example for a function of one 
variable defined at n = 0,4,8, 12, . . ., a value of the function can be estimated at a = 2,6, 10, . . . by using 
the nearest four data points. 

These techniques have been applied in a staged way to the ACSL simulation. The modified perturbation 
step size was applied to the linearization of the system, as were the unmodified perturbation step size. Then 
the data for (7 m<# was augmented by using the fourth order Lagrange interpolation for intermediate data 
points in a. Nine different a values were set, and the longitudinal eigenvalues were calculated by each of the 
three methods at each or. The differences between the various longitudinal eigenvalues were small, except 
for two cases. 'Hie locus of the eigenvalues, as a function of o was generated as a result of this work. 
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